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The relativistic corrections to the color-octet J /il) hadroproduction at the Tevatron and LHC 
are calculated up to C'(«^) in nonrelativistic QCD factorization frame. The short distance co- 
efficients are obtained by matching full QCD with NRQCD results for the subprocess g + g — 5- 
J/'4){ Sg , S\_ , Pj ) + 9- The short distance coefficient ratios of relativistic correction to LO for 
CO states ^ Sq , ^S\ ' and ^Py at large pr are approximately -5/6, -11/6, and -31/30 respectively, 
^^ ' and it is 1/6 for color singlet state ^S} . If the long distance matrix elements are estimated through 

Cn , velocity scaling rule with adopting v^ = 0.23, the LO cross sections of color-octet states are reduced 

5^ ' by about a factor of 20 40% at large pr at both Tevatron and LHC. These results indicate that rela- 

Cu ' tivistic corrections may play an important role in J/ip production and polarization at the Tevatron 

and LHC. 
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I. INTRODUCTION 



Oh. 

CD ■ 

i-i^^l Heavy quarkonium is an excellent candidate to probe QCD from the high energy to the low energy regimes. 

Nonrelativistic QCD (NRQCD) factorization formalism was established 1] to describe the production and decay of 

heavy quarkonium. In the NRQCD approach, the production and decay of heavy quarkonium is factored to short 

I> , distance coefhcients part and long distance matrix elements(LDMEs) part. The short distance coefhcients indicate 

t^^ ' the creation or annihilation of a heavy quark pair and can be calculated perturbatively with the expansions by the 

^^ strong couphng constant as- However, the LDMES, which represented the evolvement of a free heavy quark pair 

ii!^ into a bound state, can be scaled by the relative velocity between the quark and antiquark v and obtained by lattice 

^^ \ QCD or extracted from the experiment. The color-octet mechanism (COM) was introduced here. The heavy quark 

f^ ■ pair should be a color-singlet (CS) bound state at long distances, but it may be in a color-octet (CO) state at short 

^*-^ ^ distances. NRQCD had achieved great success since it was proposed. The COM was applied to cancel the infrared 

^_J . divergences in the decay widths of P-wave 0, Q and D-wave[J, Q heavy quarkonium. However, difficulties were still 

. . ' encountered. The large discrepancy between the experimental data and the theoretical calculation of J/ij] and t/i' 

^ , unpolarization and polarization production at Tevatron is a very interesting puzzle that can verify NRQCD when 

solved [GHS]. Theoretical prediction with COM contributions was introduced fitted well with the experimental data of 

«2] ! the J/ip production at the Tevatron[0|. However, the CO contributions from gluon frag mentation indicated that the 

d • J/^ was transversely polarized at large pr, inconsistent with the experimental data[lO|. 

The NLO QCD corrections of CS J/ip hadroproduction was calculated to resolve the puzzle of J/tjj 
hadroproduction JjJ, Il2l |. The calculation enhanced the cross sections at large px by approximately an order of 
magnitude. However, the large discrepancy between the CS predictions and the experimental data remained un- 
solved. The relativistic correction to CS J/^p hadroproduction were very smallflS]. At the same time, the NLO 
QCD corrections of COM J/ip hadroproduction were calculated, which eventually provided a possible solution to the 
long-standing J/ip polarization puzzle [ijj ll3 ■ A similar large discrepancy was found in the double charmonium pro- 
duction at B factories JlG'-flS'] . Much work had been performed, and apparently these discrepancies could be resolved 
by including NLO QCD corrections 19l-l22| and relativistic corrections|23-26]. The data from B factories pointed out 
that the COM LDMEs of J/tp production may be smaller than previously expected [25l - l29| . (More information can 
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be found in Refs.[30|-|32l). 

In this paper, the effect of relativistic corrections to the COM J/ijj hadroproduction at Tevatron and LHC were 
estimated based on NRQCD. The short distance coefficients were calculated up to 0{v^). Many free LDMEs were 
realized at ©(w^), which were estimated according to the velocity scaling rules within NRQCD [33|. 

The paper is organized as follows. In Section II, the frame of calculation is introduced for the relativistic correction 
of both the S- and P-wave states in NRQCD frame. Section III provides the numerical result. Finally, a brief 
summary of this work is presented. 



II. RELATIVISTIC CORRECTIONS OF CROSS SECTION IN NRQCD 

The inclusive production of J/V' was dominant by the gluon-gluon fusion process at high energy hadron colliders. 
We can compute the cross section at the level of parton. Finally, the cross section of the hadronic level is obtained 
by integrating the product of the partonic cross section and the gluon distribution function over momentum fraction. 

g(fci)+g(fc2)^J/V^(P)+g(fc3). 

As the beginning of this paper, in NRQCD, the calculation to production of charmonia is divided in two part: short 
distance coefficients and LDMEs. So the above process can be written as: 

^-^[g + g^Jl^ + g)=Y,^^{Oi'^). (1) 

at — rric 

On the right-hand side of the equation, the cross section is expanded to sensible Fock states noted by the subscript 
n. Fn, i.e., short distance coefficients, which describe the process that produces intermediate QQ in a short range 
before heavy quark and antiquark hadronization to the physical meson state. (0„ ) are the long distance matrix 

elements that represent the hadronization QQ evolutes to the CS final state by emitting soft gluons. On are local 
four fermion operators. The factor of ttiJ?""^ is introduced to make Fn dimensionless. 

In this section, our calculation on the differential cross section for this process in the NRQCD factorization formula 
is divided into two parts, namely, kinematics and short distance coefficients. 

A. Kinematics 

We denote the three relative momenta between heavy quark and antiquark as 2q^ with \q\ ^ nicV, in J/ip rest frame, 
where rric are the mass of charm quark and v is the three relative velocities between quark and antiquark. Thus, the 
momenta for the quark and antiquark are expressed as|24l.l34Ll35l| 

Pc == {Eq,q), 

p-c = {E„-q). (2) 



where Eg = \/'m?c + M^ is the rest energy of both the quark and antiquark, and 2i?g is the invariable mass of J/V"- 
When boosting to an arbitrary frame, 

Vc^\P^ q, p-c ^\P-<1- (3) 

where P is the four momenta of J/tp, and q receives the boost from (0, q). 
The Lorentz invariant Mandelstam variables are defined as 

S=(fci+fc2)2 = (P + fc3)', 
t=(fcl-P)2 = (fc2-fc3)', 
U={ki- k^f = (fc2 - Pf. 



with the relationship s + t + u = P^= AE"^. Here, when expand the ampUtude with |g| up to 0(|(f|^), great care must 
be taken to ensure that that t and u are dependent on \q\. Another two natural Lorentz invariant variables are the t' 
and u' corresponding to t and u when the effect of Iff] is considered. One can easily get that they satisfy 



s-AE^„ , 
s — Arrii 



s — Ami 



(4) 



In our subsequent calculation, we adopt t{u) to represent t'(u') directly for simplification. 

The FeynArts (36| package was used to generate Feynman diagrams and amplitudes, and the FeynCalc [37| package 
was used to handle amplitudes. The numerical phase space was integrated with Fortran. 



B. Short-distance coefficients calculation 



According to NRQCD factorization, the differential cross section of gluon-gluon fusion up to next order in v'^ to CS 
state ^S[ and each CO states ^S^ , '^S[ , '^Pj , can be expressed as 



da - ddioP^i"] + daiol'S^^^] + daiol^sf^] + daio[^pf] 



(5) 



In this expression, relativistic correction parts, denoted as ^rc\ can easily be distinguished from LO, denoted as lo\ In 
addition, each differential cross section to different Fock states should be divided in short distance coefhcient part and 
LDMEs. We can introduce F(^^+^Lj) to express the short distance coefficient of the LO cross section, corresponding 
to G{^'^~^^Lj) for relativistic correction. Many LDMEs are present, all of which are denoted by (O'^^^i^'^^^Lj)) 
and (V'^^^C^'^^^Lj)) for the LO and relativistic correction term respectively. The explicit expressions of the eight 
four-fermion operators are[l|: 



<0|C'^/^(3S'f 

< o\v'/^'epf 



|0> = <0 

|0> = <0 

|0> = <0 

|0> = <0 

|0> = <0 

|0> = <0 

|0> = <0 

|0> ^ <0 



xVV(a^av,)V'V'x|0 >, 

^[x^T^^{ala^)^^T^{~'-'^rX + h.c.]\0>, 
hx^TVtlj{ala^)i;^TV{-'-'^)\ + h.c.]\0 >, 



^-^ 



<-^ 



^fT-(--D'^anH<^la^)^^T'^i~7;D''^')x\0>, 



^-)r 



[x^T^i--'-D^an^{ala^)^^T^i-'-'^f{-'-D^a^)x + h.c.]\0> 



(6) 



Where x s-nd "0 are the Pauli spinors describing anti-charm quark creation and charm quark annihilation respectively. 
a is the Pauli matrices and D is the gauge-covariant derivative with D = D — D. T is the SU{3) color matrix. It 
should be noted that 



< 0|O^/^(3pj*l)|0 >= (2J+ 1)(1 + Oiv^)) < 0|O-^/'^'(3P(|^l)|0 



> 



(7) 



So there are 6 LDMEs for P-wave and 6 LDMEs for S-wave. The long distance matrix elements of heavy quarkonium 
decay may be determined by potential model [23, |33 or lattice calculations [40| . and by phenomenological extraction 
from experimental data[lj, |4l|. But it's very difficult for production LDMEs. To simplify the discussion of the 
numerical result, it is assumed that 



< 0\V''/^CpP)\0 > = (2J + 1) < 0|7'^/'^(3P(1*')|0 > 



(8) 



For simplicity, we can estimate the relation between their order as|38l. |3£ 



W 



(OJ/V-P^'+ii./)) 



[l + 0{v')]. 



(9) 



l^l has defined before and equal to rricV. For charmonia, v'^ ^ 0.25. Thus the total cross section can be divided in 
four sub cross section which noted by each Fock state to be calculated respectively. 

The short distance coefficients can be evaluated by matching the computations of perturbative QCD and NRQCD: 



da 



pert QCD '^m't'^^ " 



pert NRQCD 



(10) 



The covariant projection operator method should be adopted to compute the expression on the left hand of the 
equation. Using this method, spin-singlet and spin-triplet combinations of spinor bilinears in the amplitudes can be 
written in covariant form. For the spin-singlet case: 



^z;(s)fi(s)(i,s;i,s|0,0) 



{-f_ + TOc)7; 



t + 2E„ 



2^{Eq + m) 
For spin-triplet case, the expression is defined as: 



2E, 



+ rric). 



(11) 



^u(s)u(s)(-,s;-,s|l,S'^ 



-{-f_ + mc)i 



f + 2Eq 



2V2{Eq 



2E,. 



+ mc). 



(12) 



Where e denotes the polarization vector of the spin-triplet state. In our calculation, Dirac spinors are normalized as 
uu = —vv = 2mc- 

The differential cross section of each state then satisfies: 



da((2^+i)Lj)~^|X(.g + .g^(cc)((2-+i)ij)+.g)P(0^/^(2«+ii^)). 



(13) 



Where the symbol ^ is the sum of the final state color and polarization as well as the average initial states. According 
to this expression and Eq.®, expanding the cross section to 0{v'^) is to expand the amplitude squared on the right 
side of the above expression to 0(1(7]^). 

Next, we prepare to expand the short distance coefficients to the next order in v'^ . First, we expand each Fock state 
amplitude, including the S-wave and P-wave states, in term of the relative momentum \q\: 



= ^p{M? 



9=0 ^ 2'^ '^ dq'^dq^ 



q=0 



0{q'). 



(14) 



= M 



g=0 2^ ^ dq'^dqP 



q=Q 



+ 0{q'). 



(15) 



M{g + 9^{cc){^Pf)+g) = e,{s^)e„{L^){^ 



q=0 



1 . «^'(x/W^*) 



-TT?"/ 



2^ ^ dq'^dqf'dq" 



q=0 



) + Oiq'). 



(16) 



The factor . r-ir- comes from the relativistic normahzation of cc state. Odd power terms of four-momentum q vanish 

in either the S-wave or the P-wave amphtudes, where AAt and A^^ are inclusive production amplitudes to triplet and 
singlet cc respectively. 

Ms = E E(i ^' ^' ^|0' 0)(3^; 3j|l, 8a)A{g + g^c' + c^). 



,1 1 



M, =^^(-,s;-,s|l,5,)(3*;3j|l,8a)^(5 + 5^c'+c') 



In the process of expanding the amplitude squared in terms of jq], it needs to be integrated over the space angle to 
four-momentum q. We can obtain the following formulas: 



(17) 
(18) 



dfiyo*o(g)roo(g)g"/ = JkTn"^ 
J2jdnY,i^{q)Y,LA<i)q"q^ = ^kf n"^ 

T ■J 



5 

where H"'' = -g^''' + ^^. 

Subsequently, the amplitude squared up to the next order can be obtained 



^|A^(35^'1)P = A^P(0)X^(0)5]ep6l-f i|gf 



(n 



a^ _ 



dq'^dqP 



g=0 



M^*{^) 



Y^e,el + 0{v% 



Y,\MCsf)\'' = MMM:{Q) + \\cf\^ 



d\J^Ms) 



Y.\M{'pf)? = kfn 



ap 



dMt 



9g" 



dM 



A* I 



15 



kt 



(jja/3jjar 



93 



E 



fp^A 



53 



53 



.,aPtr,(y Ftr,T 



dq'^dq^dq-' dq''dq"i 



+ 0{v^). 



'-N E,'>> dqP 



(19) 



(20) 
(21) 



HE^p^a) 



+ 0{v^). 



q=0 

(22) 



Any term, which is in the order of \q\'^, must not be missed to obtain the correction up to the order of u^. In the three 
expressions above, each first term on the right side of the equation can be expressed in terms of kinematics variables 
s, t, u. Here t, u depends on \q\. The sum of terms in the order of \q\'^ in the first term as well as all the values of the 
second term is the result of the next leading order. Polarization sum ^ epe^ is equal to 11^'*'. It should be noted here 
that the polarization vector Cp also depends on q, but the dependence only appears even powers of four momentum 
q, and the explicit expressions can be found in Ref. [42]. So ep can be carried out in Eq. ipO)) . 



Therefore the differential cross section in Eq. ([5]) takes the following form 



da 
It 



ig + g^J/^ + g) 



>(35W) 



TO- 
i2 



.^O|0j/V(3^li])|O^ 



^('^0 ) (o|0J/^(i5[8])|o^ 



-{Q\0-^'^{^Sr)\Q) 



Fi'PD 



^O|0J/'A(3pM)|o) 



l + 0(f^)). 



G(35W) 



ml 

i4 



^0|pJ/^(3^[l])|0^ 



to4 
^ g('4") ^0|p.//^(3p[8]^|0A ^ 



(23) 



The explicit expressions of the short distance coefficient to the relativistic correction of CO states ^S\j and ^S^ 
are relegated to the Appendix. The short distance coefficients of CO states ^Pj are too long to list in the Appendix. 
The result of our relativistic correction of ^SJ is consistent with that of Ref.[l3| and was not given in this paper. 



III. NUMERICAL RESULT AND DISCUSSION 

We adopt the gluon distribution function CTEQ6 PDFsj43|]. And the charm quark is set as rric = 1.5Gey. The 
ratios of the CO short distance coefRcient between LO F and it's relativistic correction G at the Tevatron with 
y/s = 1.96GeV and at the LHC with -y/s = 7GeV or y^ = lAGeV are presented in Fig.l. The ratios at the Tevatron 
and at the LHC are very close. In the large pr limit, 



M2 

u 



M2 M2 

t p^ 



(24) 

Then we can expand the short distance coefficients with M . The ratios of first order in the expansion are 



i?(i4«i) 

i?(3p[«l) 



PT»M 



pt:»m 



pt:»m 



Pt^^M 



G(34il) 


F(34'') 
G(i4«l) 


G(3^fl) 


F(35fl) 
G(3p[«l) 


P(3P[«1) 



pt:>m 



pt:s>m 



pt:»m 



pt:»m 



5 
6 

11 

31 
30 



(25) 



It is consistent with the curves in Fig.l. 

The ratio of '^S[ is approximately —11/6. In large pr limit, the dominate contribution of this subprocess is 
g* — >■ cc{^S[ ). The propagator of virtual gluon g* is proportional to l/E'^. This term offer a factor of —2 to the 

ratio i?(35fl). 

As discussed in Section II, the LDMEs of relativistic correction are depressed by approximately -0.23 to LO. If 
the long distance matrix elements are estimated through velocity scaling rule with adopting w^ = 0.23, the LO cross 
sections of CO states are reduced by about a factor of 20 40% at large px at both Tevatron and LHC. Compared with 
the relativistic corrections to CS state, that ratio is approximately 4%. ^ 



1 In Ref.[l||, the ratio is about 1% for the LDMEs is different: 



(0|C'-^/'*{^5W)|0> = 1.64Gey^ (0|-P'^/'''(^ 5^)10) = 0.320GeV^ 



(26) 



The QCD corrections of both CO and CS states had been calculated in [H [Tl|. The K factor of NLO QCD 
corrections is very large for ^Pq ' and '^S\ , and it is about 1 ± 0.2 for ^S\^ and ^S\ . If we fit the matrix elements 
to NLO in as and v^ at large pt, (0|0"'/^(3S'[^1)|0) and {{)\0-'''i'{^sf^)\Q) should be enhance by a factor of 70% and 
25% respectively by relativistic corrections. But the relativistic corrections for ^P^ ' and ''^S{ can be ignored at 
large px for it is much smaller than the QCD corrections And it may affect the size of LDMEs too. 
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FIG. 1: Ratios of Oiv^) to 0{v^) results for CO states. The left diagram's result is calculated at Tevatron and the right one 
is at LHC. They give pretty close values at large pr- 



IV. SUMMARY 



In summary, we calculate the relativistic correction terms to CO states for J /^ production at the Tevatron and at 



the LHC. The short distance coefficient ratios of relativistic correction to LO for CO states ^5, 



[8] 3c.[8] 



^S^P and 3p 



at 



large pt are approximately -5/6, -11/6, and -31/30 respectively, and it is 1/6 for color singlet state ^S{ . If the long 
distance matrix elements are estimated through velocity scaling rule with adopting v^ — 0.23, the LO cross sections 
of CO states are reduced by about a factor of 20 40% at large pr at both Tevatron and LHC. Compared with the 
relativistic corrections to CS state, that ratio is approximately 4%. Thus the result may affect the production and 
polarization of J/ij} and hadronic colliders. And it may affect the size of LDMEs to. 

These results indicate that relativistic corrections may play an important role in J /if) production and polarization 
at the Tevatron and LHC. 



Then 



which is much smaller than v^ 



(i:)\V'/^{^S^^^)\0)/{0\O-'/'^{^S^^^)\Q)/ml = 0.087 



(27) 



■ 0.23. 
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V. APPENDIX: SHORT DISTANCE COEFFICIENTS 



The short distance coefficients of ^5n were: 



ml 



320 



1111 



167rs2 64 4 Ar2 - 1 
M^2 (t^ + tu + u^) + [i^ + tu + u^f - M^° (4i^ + ifu + Itu^ + 4u^) 



+M^ {8f + SSt'^u + 66t\^ + Slt^u^ + mfu^ + 33iu^ + 8m^) 



M'^{f + tu + u^) {4£^ + %fu + %tu^ + Au^) 



M {M^ - tft (M^ - uf (M^ -t-u)u{t + uf 



(28) 



G{'sf^) 



1 11 
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m^ 



640 



167rs2 64 4 iV2 - 1 
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280tu^ + 95u'') 



+M^^ (95t'' + 280t^u + 358t^i 

-3M*2 (l6i^ + 9m^u + 175t^u^ + 175i^u^ + OStu"* + 16u^) 

-2M^° (45i*^ + 72f^M + 21t^u^ - 22t^u^ + 21t^u^ + 12tu^ + 45^*^) 

+M® (I98i^ + 678t®u + 1141*^^^ + 1345^*^^ + 1345t^M'* + 1141*^^5 + 678iM^ + 198u^) 

-M^^80t^+756t^M+1583t^uV2224t''^M3+2446t'*M'* + 2224f\'^ + 1583^^^^ + 756tu^ + 180m® 

+M'*(85i^ + 408t*^M + lOOOt^M^ + 1637i^7? + 2028i^/ + 2028t'*7i''^ + 1637t^M*^ + lOOOt^w^ 

+408tM*^ + 85w^) - M^ {t^ + 2t^u + 2tu^ + u^Y {llf^ + 30f''M + 30^^^^ + 30tM-'' + 17?/) 
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The short distance coefficients of '^S\ were: 
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(30) 



3Q[8]^ 



G(3^1' 



1 1 1 



1 



167rs2 64 4 7V2 _ 1 3 



(-512) 



Af " (87t2 + 22tu + 87^2) + M^^ (_i4^3 ^ 335^2^ ^ SSStu^ - 14u^) 



2^//'*^ 399r + 1612ru + 2020ru^ + 1612fu'' + 399u 



2„,2 



,3 



-M^ (2100^^ 



8976ru 



14497t3u2 j^ I4497t2u3 . 

M^ (2590t^ + 12096t^u + 23855*^^^ + 2Q?,lU^u^ + 23855t^u 
^- ' 19905tV + 29152tV + 29152iV 



8976to* + 2100u^ 

2„,4 



+M^ (l620i^ + 8498t-u 



^27M2 (l5t^ 



lOifu 



295t^u^ 
2 , J... , ,.2 



+297iM(i + u) (r + iu + -u' 




612ru 



4„,4 



- 12096iu^ 
19905t2u5 - 
blOt^u^ + 295fu^ 



tf (Af2 



ufit- 



■u 



- 2590u''} 

8498ta^ + 1620u^) 
^ 104to^ + 15m^) 



(31) 
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